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The observational consequences of certain broken symmetries in a thermalised quark gluon plasma 
are elucidated. The signature under study is the spectrum of dileptons radiating from the plasma, 
through gluon fusion. Being a pure medium effect, this channel is non- vanishing only in plasmas with 
explicitly broken charge conjugation invariance. The emission rates are also sensitive to rotational 
invariance through the constraints imposed by Yang's theorem. This theorem is interpreted in the 
medium via the destructive interference between various multiple scattering diagrams obtained in 
the spectator picture. Rates from the fusion process are presented in comparison with those from 
the Born term. 
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I. INTRODUCTION 



fSJ ' Experiments are now underway at the Rclativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 

^ ] to study nuclear collisions at very high energies. The aim is to create energy densities high enough for the production 

OO • of a state of essentially deconfined quarks and gluons: the quark gluon plasma (QGP). The QGP is rather short lived 

[ and soon hadronizes into a plethora of mesons and baryons. Hence, the existence of such a state in the history of a 

^ . given collision must likely be surmised through a variety of indirect probes. One of the most promising signatures 

I ■ has been that of the electromagnetic probes: the spectrum of lepton pairs and real photons emanating from a given 

, collision. These particles once produced interact only electromagnetically with the plasma. As a result they escape 

' the plasma with almost no further rescattering and convey information from all time sectors of the collision. 

In this article, the focus will be on the spectrum of dileptons radiating from a heavy-ion collision. The primary 

^ • motivation for measuring such a spectrum is the hope that the formation of a QGP in the history of a collision 

I ' will produce a qualitative or quantitative difference in the observed rates. The measured quantity is the number of 
dileptons, usually binned according to their invariant mass. It is assumed that this may be estimated by the following 

^ ' factorized form, 

> : 

■ dN r^f r^+it) r^+it) rv+it) r M d^E 

X ■ = / / / / / d\^^''-^^{q\lnt,xl^J.{t,x)). (1) 

■ dM J^^ J q^> d-^q 

where, — Jiq"' the number of lepton pairs produced per unit space time, per unit four-momentum, from the unit 
cell at {x,t) in a plasma in local equilibrium (local equilibrium is assumed here). Ostensibly, this depends on the 
four-momentum of the virtual photon {q^,q), the temperature (T), and the relevant chemical potential (/x). The 
temperature and relevant chemical potential are, in general, local properties for an expanding plasma and vary from 
point to point in the plasma as indicated. Finally, the rates from each space time cell have to be integrated over the 
entire space time evolution of the plasma; where, the spatial limits of the expanding plasma are represented by the 
variables x-{t),x+{t), y-{t),y+{t), z-{t),z+{t). 

Many calculations of the dilepton radiation in the deconfined sector have concentrated on the Born term qq e~^e~ 

to estimate the differential rate ^S"^"" Q • In those, the focus has usually been more on the effect of the space time 
evolution of the plasma on the final spectrum. Higher order rates have also become recently available 01- AH these 
rates essentially consist of vacuum processes that have been generalized to include medium effects of incoming medium 
particles along with Pauli blocking (Bose enhancement) for outgoing fermions (bosons). They also include thermally 
generated widths and masses for the propagating particles. However, these rates have non-zero vacuum counterparts. 
Contrary to these are the pure medium reactions, i.e. processes whose vacuum counterparts are identically zero. Such 
processes arise as a result of the medium breaking various symmetries which are manifest in the vacuum jl^ . The 
motivation behind exploring such processes is the possibility of observing a spectrum (emanating from these) which 
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FIG. 1: Diagrams that are zero by Furry's theorem and extensions thereof at finite temperature. These become non-zero at 
finite charge density. 



is noticeably distinct experimentally from the case where a QGP was not produced in a collision, or the symmetry 
remained unbroken by the plasma. Such a channel will be explored in this article. 

It is now well established that the central region at RHIC is not just heated vacuum, but actually displays a finite 
baryon density Q or an asymmetry between baryon and anti-baryon populations. This asymmetry may be achieved 
by the introduction of a quark (or baryon) chemical potential fiq. For example, it may be argued that any baryon 
number asymmetry prevalent in the QGP must have been introduced by valence quarks, which, having encountered 
a hard scattering, failed to exit the central region. Hence, a chemical potential fi is provided for the up and the down 
quark. The strange quarks are brought in by the sea or produced thermally in the medium in equal proportion with 
anti-strange quarks. Hence, they are assigned /i = 0. In most heavy-ion collisions, the nuclei of choice are rather large 
and display isospin asymmetry, i.e. there is an asymmetry in the populations of neutrons and protons being brought 
into the central region. If the stranded valence quarks in the plasma arrive with equal probability from either nucleon, 
one would require a higher fi for down quarks. As a first approximation, this effect is ignored, and, in the remaining, 
accept fXu = fJ-d- 

In such a scenario a finite baryon density may lead to a finite charge density, this is discussed briefly in Sec. II. It 
has been proposed that the presence of a finite charge density will lead to a new channel for the production of lepton 
pairs 01 . Diagrammatically, this is achieved through a two-gluon-photon vertex with a quark triangle (see Fig. [5J). 
The vacuum counterpart of this process is constrained by Furry's theorem and is identically zero. The extension of 
this symmetry and its breaking by the medium was discussed in 0; for completeness a discussion is included in Sec. 
II. Calculations here will be carried out in the imaginary time formalism however our method of treating finite 
density will differ slightly from the standard method. This is outlined in Sec. III. Our formalism leads naturally to the 
spectator interpretation Q of the quark loop, also discussed in Sec. III. In vacuum, gluon fusion is also constrained 
by Yang's theorem [lo| . This constraint, based on rotational invariance also sets the vacuum counterpart to zero. 
The extension and eventual breaking of this symmetry in the medium are discussed in Sec. IV. Yang's theorem is 
broken by two different medium effects: each is isolated and the dilepton rate from it is evaluated in Sees. V and VI. 
Concluding discussions are presented in Sec. VII. A brief appendix of intermediate derivations follow. 

II. BARYON DENSITY, CHARGE DENSITY AND FURRY'S THEOREM 

At zero temperature, and at finite temperature with zero charge density, diagrams in QED that contain a fermion 
loop with an odd number of photon vertices {e.g., Fig. ^ are canceled by an equal and opposite contribution coming 
from the same diagram with fermion lines running in the opposite direction, this is the basic content of Furry's 
theorem (see also This statement can also be generalized to QCD for processes with two gluons and an 

odd number of photon vertices. The theorem is based solely on charge conjugation invariance of the theory. 

In the language of operators, it may be noted that these diagrams are encountered in the perturbative evaluation 
of Green's functions with an odd number of gauge field operators i. e. 



(0|^plA^2-.^M2n+l|0). 
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In QED, CA^C ^ = — A^, where C is the charge conjugation operator. In the case of the vacuum, C|0) = |0). As a 
result 

= (0|A^,A,,...A,,„,jO)(-l)2"+i 

= -m^,A,,,...A,,,^^,\Q)=Q. (2) 

In an equihbrated medium at a temperature T, we not only have the expectation of the operator on the ground state 
but on all possible matter states weighted by a Boltzmann factor i.e. 



where (3 = 1/T and ^ is a chemical potential. Here, C\n) = e*"^] — n), where | — n) is a state in the ensemble with the 
same number of antiparticles as there are particles in \n) and vice- versa. If /Lt = one obtains 

{n\A^,A^,...A^„^^,\n)e-^^" = -{-n\A^,A^,...A^,^^,\ ~ n)e-^^'\ (3) 
The sum over all states will contain the mirror term (— — n)e^^^", with the same thermal weight 

5^(n|A^,A^,...A^,„^Jn)e-^^" = 0, (4) 

n 

the expectation over states which are excitations of the vacuum |0) will again be zero as in Eq. and Furry 's 
theorem still holds. However, if /i 7^ 

(n|A^,A^,...^^,„^Jn)e-^(^"-^'3") = -{-n\A,,A,,...A,,^^,\ - n)e-^^^--^^-\ (5) 
the mirror term this time is (— ri| ...A^^^^ J — n)e^^'^^^"^^'^"\ with a different thermal weight, thus 

^(n|A^,A^,...A^,„,Jn)e-^(^"-^«") ^ 0. (6) 

n 

This represents the breaking of Furry's theorem by a medium with non-zero charge density or chemical potential. 

Some points are in order: there is more than one kind of density that may manifest itself in the plasma. There 
is the net baryon density which requires that there be a difference in the populations of quarks and antiquarks of a 
given flavour. There is the net charge density which simply requires that there be more of one kind (either positive 
or negative) of charge carrier in the medium. Note that it is possible to have a net baryon density and yet no charge 
density and vice- versa as Tabled indicates. As mentioned in the introduction, it will be assumed that there is a net 
baryon density, which manifests itself solely in the up and down flavours of the quarks. As the up quark has a charge 
of -f| and the down quark — i; equal densities of both will lead to a QGP with a net electric charge density. It will 
be demonstrated in the following sections that it is this density that leads to a dilepton signature of the breaking 
of Furry's theorem at leading order in the EM coupling constant. The baryon density merely serves the purpose of 
generating such a charge density. Hence, this signal is not present, at leading order, in a plasma with a = /j.^ = yUg, 
where the net charge is zero [23 . 

In the previous paragraph, pure QED diagrams have been discussed . One may now make the most simple extension 
to QCD, by replacing two of the photons with incoming gluons. It is to be noted that while the photon is an eigenstate 
of the charge conjugation operator C the gluon is not There are eight gluons, each carrying a colour charge in 
the adjoint representation of SU(3). The sole role played by colour in this calculation will be to furnish the factor of 
Tr[f°t''] in the Feynman rules. The calculations arc identical to those in QED. The reasons for considering this sort 
of diagram over others are obvious: this is the lowest order effect in the series, loops with more particles attached 
will invariably be suppressed by coupling constants and phase space arguments. Also, diagrams with more than two 
gluons are non-zero in the vacuum itself and finite density effects may then be a mere excess on top of an already 
non-zero contribution. Our exploratory calculation mainly seeks to highlight the behaviour of a new channel. 

Cases with all values of the three-momentum p of the photon from zero (maximum timelike) up to almost the 
energy E of the photon (almost lightlike) will be considered. We will consider cases where the gluons will be both 
massive and massless. The quarks will be massive in all cases. 
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TABLE I: Different scenarios of plasmas with different baryon and charge densities. 
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FIG. 2: The two gluon photon effective vertex as the sum of two diagrams with quark number running in opposite directions. 



III. FORMAL CALCULATION 

In this section, the computation of the dilepton production rate from the two-gluon channel is outUned. The first 
step is the evaluation of the two-giuon-photon vertex in the imaginary time formalism. This is then used to construct 
a three-loop photon self-energy. The imaginary photon frequency is then analytically continued to real values. On 
the real axis one encounters various branch cuts: we pick the cut that corresponds to the process of two-gluon fusion 
and evaluate the imaginary part of the photon self-energy. This is a rather technical procedure. However, a method 
has been proposed which allows one to construct the imaginary part of the photon self-energy in terms of multiple 
scattering diagrams This technique has been loosely termed as the "spectator interpretation" of self-energies. A 
detailed investigation of this procedure has been carried out for two loop self-energies in cj)'^ theory as well as in 
QCD 01 • In the following subsections the spectator interpretation will be extended to three loops. This represents 
a significant application of the spectator interpretation, and this permits a physical explanation of the extension and 
breaking of Yang's theorem in the medium. 

A. The two-gluon-photon vertex at 7^ and p 7^ 

In the following, we will outline the formal derivation of the two-gluon-photon vertex. Details of the method are 
presented in the appendix. The Feynman diagrams for the two-gluon-photon vertex as illustrated in Fig. [21 consists 
of two sets of quark triangle diagrams with the fermion number running in opposite directions. The two vertices are: 
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}i + m) 



q2 — 



(q-k)2 



(q — p)2 — 



(7) 



3 J {2Tif 2^ V ^ '''^ (q + p)2-m2' 



(q + k)2 — m? 



q2 — 



(8) 



where the trace is imphcd over both colour and spin indices. As always in the imaginary time formalism, the zeroth 
components of each four-momentum is a discrete even or odd frequency: 



i{2n+l)TTT + 



P 



i2m7rr, 



i2jTTT. 



n, m, j are integers, n is the quark chemical potential, and the are Gell-Mann matrices. The overall minus sign 
is due to the fermion loop. The sum over n runs over all integers from — oo to +oo. This sum may be performed 
by two distinct methods: the method of contour integration [g and the method of non-covariant propagators [l5l |. 
Each method is more advantageous in certain cases. In this article, we use the method of contour integration to 
evaluate Eqs. {T)) and (jSJ (for an evaluation of similar diagrams using non-covariant propagators see Ref. Q)- For 
later convenience, we will separate the momentum dependent and mass dependent parts of the numerators of both 
j-fii^p Q^jj^ T'MP'^ ^ This is merely a formal procedure and for T^^f consists of the following: 
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2(3 J (27r)3 
(q2 — m2)((q — k)^ — m?){{(\ — p)^ — i 



-I- m 



.A^'^^Pqa + AP-^'I^Piq - k),3 + A^"'P'^iq - p)- 

(q2 — TO^^^^q _ Jj^2 _ 77i2)((q _ p)2 _ 



(9) 



Where A^''^^'^ represents the trace of four 7 matrices and 3^"°"^^^^ represents the trace of six 7 matrices. The 
denominators of both Tia^p-y and Ai'^'' are the same and hence have the same set of poles. 

The sum over n may be formally rewritten as a contour integration over an infinite set of contours each encircling 
the points q^ = i{2n + l)7rT -I- /i. The difference between this situation and that at zero density (/i = 0) is that the 
contours are on a line displaced by from the y axis. In the usual procedure (see Sec. (3.6) of Ref. Q), one separates 
the vacuum piece, a thermal particle and antiparticle piece, and a pure density contribution. However it is possible to 
deform the contours in a way entirely similar to the zero density situation. One obtains two infinite closed semi-circular 
contours: see appendix for details. The sole difference from the zero density situation is that one of the contours will 
be multiply connected: in the case of T^^p this consists of the exclusion of the points at g° = i{2n + \)'kT + fihy 
an infinite set of infinitesimal contours, while in the case of T^p'^ ^ the points at q^ = i{2n + VjiiT — fj, are excluded. 
Performing the contour integration will essentially result in the evaluation of the residues of the functions in Eqs. Q 
and (jS)) at its various poles with appropriate finite density thermal weights. Combining the results obtained from the 
application of this procedure on T^-'^p and T^'"' one obtains 



2/3 



-Res. 



6/3(9°-^) + 1 e'3(«"+A') + 1 



(q2 — m^){{c[ — k)2 — ?T!,2)((q — p)^ 



4m- 



.gt'^iq -p-k)P + gt'Piq - k + p)" + g''P{q + k - pY 
(q2 — w?){{(\ — k)2 — m?)[[<\ — p)2 — m?) 



1 



1 



(10) 
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Where, the w^'s (with i running from 1 to 6) are the residues of the function within the large square brackets. We find, 
as would have been expected, that the entire contribution is proportional to the difference of the quark and antiquark 

distribution functions. We denote these as An(q^,ri) ~ ( „,,)"'" , „, oi . — ). The residues will be evaluated at 

the various poles of the integrand. A close inspection of Eq. (|10|l indicates that there are three poles on the positive 
X axis at, 



^/fT^ = Eq (11) 



= J\q-k\^+m^ + k" =Eg-k + k'' (12) 



= ^\^-P\2 + m^+p^ ^Eg^,+p'. (13) 



and three on the negative x axis. 



= -^q^+7n^ = -Eg (14) 



(?" = -J\q-k\^+m^ + k° = -E,-k + k° (15) 



-./W^W+^ + P° = -E,-v+p''. (16) 



We denote the residue at each of these poles as residues (1-6). Before evaluating the function at each of these residues, 
we consider the fate of the remaining imaginary frequencies in the expressions k^ . The even frequency k^ also has 
to be summed in similar fashion as . The external photon frequency will have to be analytically continued to a 
general complex value and finally the discontinuity of the full self-energy across the real axis oip^ will be considered. 
We perform this procedure in the next section. 



B. The photon self-energy and its imaginary part 

We are now in a position to calculate the contribution made by the diagram of Fig. |21 to the dilepton spectrum 
emanating from a quark gluon plasma. To achieve this, we choose to calculate the discontinuity of the photon self- 
energy as represented by the diagram of Fig. |21 across the real axis of p'^ . In the previous section we wrote down 
expressions for T>^'^''(p, k,p — k): the vertex with the the two gluon momenta incoming and the photon momentum 
outgoing. To write down the expression for the full self-energy we also need expressions for T"^'^''{—p,—k,k — p): 
the vertex with the photon momentum incoming and the gluon momenta outgoing. This vertex also admits a 
decomposition into two pieces for quark number running in opposite directions. 



T'^^'i-p, -k, k-p)^ T^""' + T'^"" 



where the factor T'^^'' can be written as 



(17) 



2xbc 



q 



yP" ■ 2(3 J (27r)3 

Bt^'y-Pp'-q^{q + k)p{q+p)^ 
(q2 — 7Ti2)((q + k)2 — m?){{(\ + p)^ — m? 

A'^'^^Pq-y + A'^'^Piq + k)p + A^'''P"{q + p)a 
The traces of four and six 7 matrices admit the following identities: 



E 



(18) 



j^P-vPp 
j^p.upa 



j^pu-tp. _ j^ppu-t 
j^pPvp _ j^ppliv 
j^apvp _ j^papv 
j^apl3i^yp j^papPuj 



(19) 
(20) 
(21) 
(22) 
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In each equation above the first equaHty uses the fact that the trace of n 7 matrices in a particular order is the same 
if the order is fully reversed. The second equality uses the cyclic properties of the trace to put 7'' at the start in each 
case. Substituting the above identities in Eqs. (fH^ and ifTTji . wc may easily demonstrate that 

T'^^'i-p, -fc, k-p)^ T^-'Pip, k,p~k) (23) 

Implementing the above simplifications we may, formally, write down the full expression for the three-loop photon 
self-energy in the imaginary time formalism as (note, there are many other photon self-energy diagrams at three loops, 
however we only consider the contribution which arises due to the breaking of C invariance) 




FIG. 3: The Full photon self-energy at three loops and the cut that is evaluated. 



The diagram that we are considering is that of Fig. |3| In the above equation Vp^ik) is the gluon propagator. We 
perform the calculation in the Feynman gauge for the gluons: 

-Dfdk) = (25) 

In order to calculate the differential rate of dilepton production we need to evaluate the discontinuity of the photon 
sclf-encrgy. This involves, first, converting the sum over discrete fc*^ frequencies into a contour integral over a complex 
continuous fc", as was done for g". This would be followed by the evaluation of the contour integral by summing 
over the residues of the integrand at each of the poles of fc". Finally we look for poles and branch cuts in this 
expression as a function of p^ by analytically continuing p'^ onto the real axis. There are many poles in fc*^ for which 
we have to evaluate residues. Some of these poles are in the denominators of the gluon propagators, while some 
are in the vertices T^'^p and T'^'^^. As the residues at each of these poles is analytically continued in from a 
discrete imaginary frequency to a complex number and finally to a real continuous energy, various branch cuts will 
appear. This procedure of evaluation of residues and analytic continuation may also be performed prior to the (Pk 
or (Pq integrations: at this stage the branch cuts on the realp*^ axis appear as poles in the d^k or (Pq integrations as 
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p —^ E + ie. The presence of the ie will allow each integrand to be unambiguously broken up into a set of principle 
values and imaginary parts. Combining these will lead to various real and imaginary parts of the full self-energy and 
will correspond to various physical processes of photon propagation and decay in the medium. Twice the integral 
over the imaginary parts will give us the required discontinuity. 

To obtain the discontinuity we essentially choose a pair of poles in the expression; evaluate the residue of the dfc" 
integration at the first pole and twice the imaginary part of the cPk,d^q integration as ^ E + it at the second 
pole. Each such combination constitutes a 'cut' of the self energy or a part of a cut of the self-energy. The cut line 
essentially passes through a set of propagators in the self-energy dividing it into two disjoint pieces (see, for example, 
Fig. O. The propagators that have been cut are indicated by the energy momentum delta functions obtained from 
the residue and discontinuity procedure. If we denote the Feynman rule for one of the disjoint pieces as A^i and the 
other by M.2 then this particular discontinuity of the self-energy gives the Feynman rule for A^2-^i oi' M.iJ^2- If 
the cut is symmetric i.e. Mi = M2, then we obtain the square of the amplitude for the process |A^ip. For this 
calculation we are solely interested in the square of the amplitude of the process shown in the lower panel of Fig. |31 
Our preceding discussion indicates that this will be given by the cut line indicated in the upper panel of the figure. 
This is the process of gluon giuon fusion to produce a heavy photon resulting in a dilepton. The other cuts represent 
extra finite density contributions to processes already non-vanishing at zero density. 

The above discussion indicates that we merely have to look for poles in the denominators of the gluon propagators. 
Isolating this piece form Eq. H24|l . we note the denominators of the two sets of gluon propagators is 



k2(p-k)2 - k)ik° + k) ip" - - Ep-k)ip° - k° + Ep-k) ' ' 

Where Ep^k = \p — k\- The integration will encounter four possible poles at fc" = ±fc, and = p° ± Ep-k- Each 
choice will lead to a different process as is analytically continued to the real axis. All choices will not lead to the 
desired process. We now investigate each of these possibilities in turn. 

We begin by evaluating the residue of the remaining integrand at the pole k'^ = k. At this pole the remaining three 
denominators are 



1 1 1 

2k pO -k^ Ep_k p°-k + Ep_k ■ 

On analytically continuing p^ we will obtain two locations on the real line of where a discontinuity may occur: 
p^ = E = k + Ep-k and p^ = k — Ep-k- The second pole will lead to the photon invariant mass E"^ — p^ < Q i.e. a 
spacelike photon, we ignore this cut. Substituting the first value for p^ , wc obtain the discontinuity of the self energy 
a± E ^ k + Ep^k- This turns the gluon denominators into 



-inSiE - fc - Ep^k)-} 



2k 2Ep-k 

Evaluating the residue at fc*^ = — fc, we obtain the remaining denominators as 



-2k -39O -k- Ep^k -k + Ep_k 

This leads to two possible locations on the real line of where a discontinuity may occur: — fc — Ep-k and Ep-k — k. 
The first choice leads to a negative energy photon and the second to a photon with a spacelike invariant mass, thus 
we ignore this fc° pole altogether. 

Evaluating the residue at fc*' = p*^ + -Bp-fc and analytically continuing we once again obtain a negative energy 
photon and a spacelike photon and thus this residue is ignored as well. The final residue is at fc" = p*^ — Ep^k- This 
leads to possible discontinuities at p'^ — E = k + Ep^k and Ep-^ — k. The second possibility leads to a spacclikc photon 
and is ignored. The first gives a timelike photon with positive energy and thus is included in the cuts considered. 
With this choice we obtain the gluon denominators into 

-iTT5{E - k - Ep^k)-} 



2k -2E, 



Thus, in performing the sum over the Matsubara frequencies fc" we will only confine ourselves to two poles: one on 
the positive side of the real axis at = k, one on the negative side at k^ = p^ — Ep^k- For both poles, we analytically 
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continue p'^ to E = k + Ep^k, leading to 



p" — Ep_t^ — E - Ep_k — k + Ep_k - Ep_ 



k 



Thus, in the rest of the expression we wiU simply replace k^ ^ k and use the appropriate distribution functions in 
each case depending on whether the initial pole was on the positive or negative side. Then we will use the delta 
function to set the value of k. The results of this procedure as well as the final expressions and their properties will 
be discussed in the next subsection. 

One may also expect the gluons to acquire a thermal dispersion relation in the hot QCD medium (sec Ref. 
In a later section we will employ a simplified version of the in-medium gluon dispersion relations: the gluons will be 
ascribed a thermal mass. The above derivation of the photon self-energy and the pole analysis are still valid provided 
we use a massive vector propagator such as: 



9pC 



^PCW=-^^. (27) 

and we substitute in the vertex expressions every occurrence of the massless gluon energy k by its massive equivalent 
Ek = \ k"^ + TO? where mg is the thermal gluon mass. 



C. The spectator interpretation 

In this section we evaluate the particular cut of Fig. of the three-loop photon self-energy. Focusing on the two 
poles of k^ highlighted in the preceding subsection and performing the associated analytic continuation of we obtain 
the discontinuity in the photon self-energy as 



Discn^"^ = 



^T^'-^iE, k, Ep^,f-^T<\E, k, Ep.,)r^ 
'I 



1 



2 e/^*; - 1 
d^k 



2k ^ ^ ^ '^'2Ep^k 

{-l)(^-27:iS{E-k-Ep^k)) 

i-l)(-27:iS{E-k-Ep.k) 



1 



1 



Combining the gluon distribution functions and using the relation Ep^k + k = E, wc obtain 



(28) 



DiscH'''' 



-f^Tt^PVE, k, E„^k)^T<\E, k, E^_k)-^ 



X e 



1 ]n{k)n{Ep.k){2m5{E -k- Ep_k) 



p—k 



(29) 



Disc[-m'^''l. We 



To obtain the differential rate for dilepton production we need the quantity r = — .9^1/ 

substitute the expression for Disc [11^"] and note that the intervening factors of the metric as well the factor g^u 
are obtained from the sum over the polarizations of the incoming gluons and the outgoing photon: 



^£,%(k)£,ak)-- 



5pC 



(30) 



(31) 
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Substituting the above relations into r, we obtain 



EEE 

I j I 



<pk 



(2TTf2k2Ep-k 



X (^e'^^ - 1^ n{k)n{Ep^k) (2Tr6{E -k- Ep^k)) 



(32) 



Introducing factors of 2t: and extra delta functions we may formally write the above as a straightforward kinetic 
theory equation: 



X (^e^^ - 1^ nik)niEp-k) {{27:)^S\p - k - w)) , (33) 

where Mij.i is the thermal matrix element for two gluons in polarization states i,j to make a transition into a photon 
in a polarization state The entire process is weighted by the appropriate thermal gluon distribution functions and 
has the usual energy momentum conserving delta function. In this calculation both gluons are massless; thus they 
have only two physical polarizations. The photon being massive has an extra polarization Eg. Note that the thermal 
matrix elements Mij.i still contain thermal distribution functions, they thus encode information regarding incoming 
and outgoing quarks from the process into the medium. Using the polarization vectors, the thermal matrix elements 
may be easily expressed as vacuum multiple scattering diagrams with thermal weights for the incoming and outgoing 
quarks as well. This reformulation of the thermal loops is the spectator interpretation. We now no longer have a 
quark loop: it is replaced with a set of coherent tree diagrams. 

Recall that in the evaluation of T'^'^p we had performed a contour integration over q'~' and obtained six residues (see 
Eqs. (|ll|) - (|16|l ). We did not elucidate the residues at the time, as wc still had two complex frequencies, and p", in 
the expressions. Once the residue in the fc*^ integration has been taken and p'^ analytically continued, we obtain the 
following results. 

For the pole at = ^/q^ + rn^ = Eg we obtain. 



2rfcc 



(27r)3 

2£;g((q-k)2 -m2)((q-p)2 -m2) 
.g'^'iq -p-k)P+ gt'Piq -k + pf + g^Piq + k - pY 



2£;,((q-k)2-m2)((q-p)2 



^n{E,,ii) (34) 



Once again B'^J' represents the trace of six gamma matrices. The term above may be reinterpreted as 



d^q 



{2iTf2Ec 



'^1 



1 - h{Eq,y) - n{Eg,^) 
'" '"l - n{Eq,p) - n{Eq,n) 



((q — k)2 — m'^){{q — p)^ — 
1 - h{Eg,-p.) - n{Eq, -fi) 

1 - h{Eq, -p) - n{Eq, -fi) 



(35) 



Where r is the spin of the quark (or antiquark) of momentum q. The distribution functions have been written in 
a way to distinguish the contributions from quarks and antiquarks. If we concentrate only on the coefficient of the 
quark part of the distribution functions wc note that fi] ^ ^ = Si ^t'^''^ e* ^e* ^ is simply the Feynman rule for the process 



I p J- 



indicated as Tl in Fig. ^ The spins for the incoming quark have been averaged over, while its momentum has been 
integrated over. One may also show that the coefficient of the antiquark part of the distribution functions corresponds 
to the diagram referred to as T4 in Fig. ^ with the in coming quark line replaced by an incoming antiquark line. 
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FIG. 4: Multiple scattering expansion of the quark loop. Each diagram corresponds to a residue of the q° integration. 



Following the procedure, as outlined above, one may easily demonstrate that each residue of q'^ corresponds to a 
multiple scattering topology. As a result there are six different coherent tree diagrams as shown in Fig. 01 Each tree 
diagram in Fig. ^corresponds to a residue of the q'^ integration. No particular time ordering is implied except that 
the gluons arc incoming and the photon is outgoing. The quarks can be both incoming and out going. As a result 
the thermal matrix element Mi.j^i in Eq. H33|l may be expanded as 



{1 - n{Eq,i^L)} - n{Eq,fl) 
{1 ~ n{Eq, -/i)} - h{Eg, 



(36) 



Where Miji is the vacuum amplitude of the diagram referred to as Tl in Fig. ^ and fl} j i is the same diagram 
with the in and outgoing quark replaced with an antiquark. It may be easily demonstrated, by a simple variable 
transformation, that flj j ; = —fJ-jj i- The same is true for the other amplitudes with antiquarks, each is the negative 
of a different vacuum amplitude from the six. As a result we obtain: 



M 



4 



-h{Eq,fi) +n{Eg,-ii) 



(37) 



From the above expression it is obvious that if the chemical potential /i = then the rate is zero as well. Note 
that the imaginary time formalism only provides us with the square of the above term. Indeed it is A^^ which will 
eventually determine the rate. The uncoupling of the to individual tree amplitudes constitutes the proof of the 
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FIG. 5: Multiple scattering expansion of the lower diagram in Fig. El Each dia gram corresponds to a residue of the q 
integration. 



spectator interpretation of the imaginary time formalism. One may derive the rate starting from these simple tree 
amplitudes, without invoking the complicated machinery of the imaginary time formalism. 

We may now state the spectator interpretation for the square of matrix element. This is shown in Fig. [S] This 
is the spectator interpretation of the loop diagram of Fig. O It represents the process of two gluons in states (i, j) 
encountering two incoming medium quarks (or antiquarks) with quantum numbers q, Q leading to the emission of 
a photon in state / and two quarks (or antiquarks) with identical quantum numbers q, Q. In the amplitude on the 
left hand side of Fig. |31 q participates in the reaction whereas Q is a spectator. In the amplitude on the right the 
reverse is true. Note that we do not require q, Q to be simultaneously quarks or antiquarks, they may be either. We 
have thus expressed the complicated loop containing matrix element as a coherent sum of simpler tree diagrams. The 
main purpose of such a decomposition is more than just a physical perspective: it allows us an understanding of the 
mechanism of symmetry breaking not provided by the rules of the imaginary time formalism. This will be discussed 
in the subsequent section. 

In passing, we should once again point out how the spectator interpretation greatly simplifies any thermal cal- 
culation. The diagrams of Fig. 21 are not difficult to motivate from first principles. They represent the set of all 
possible means (at lowest order in coupling) by which one may couple two gluons to a photon with two fermions from 
the medium. Along with this is the restriction that the fermions return back to the states that they vacated. The 
six diagrams represents the six different means of ordering the gluons and the photon. The two kinds of thermal 
factors {h(Eq,±fi), {1 — n(Eq,±fj,)}) represent the possibilities of the fermion being ejected from the medium prior 
to its re-absorption, and vice-versa. The two fermions are allowed to be both quarks and antiquarks as required by a 
relativistic medium. 
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IV. ROTATIONAL INVARIANCE AND YANG'S THEOREM. 



The vacuum analogue of the two-gluon-virtual photon process does not exist due to Furry's theorem. If it did, it 
would represent an instance of two identical massless vectors fusing to form a massive spin one object; or alternatively 
a massive spin one object decaying into two massless vectors. There exist other such processes not protected by Furry's 
theorem, e.g., to — > 77. Such a process though not blocked by Furry's theorem is still vanishing in the vacuum. We 
effectively have a situation where there are two massless spin one particles in the in state and a spin one particle in 
the out state, or vice versa. In such circumstances another symmetry principle is invoked. This symmetry principle, 
due to C. N. Yang is based on the parity and rotational symmetries of the in and out states and will, henceforth, 
be referred to as Yang's theorem. 



A. Yang's theorem in vacuum 



The basic statement of Yang's theorem, as far as it relates to this calculation, is that it is impossible for a spin one 
particle in vacuum to decay into two massless vector particles. This statement is obviously also true for the reverse 
process of two massless vectors fusing to produce a spin one object and as a result a fermion and an antifcrmion 
combination in the triplet state. This may be understood through the following simple observation. Imagine that we 
boost to the frame where the two incoming vectors (in this case gluons) are exactly back-to-back with their three- 
momenta equal and opposite. The outgoing vector (the virtual photon in this case) is produced at rest and eventually 
disintegrates into a lepton pair. We will now apply various symmetry operators (parity, rotation, etc.) on both the 
incoming and outgoing states. Note that, as we arc only interested in strong and electromagnetic interaction, hence, 
parity is a good quantum number. If both incoming and outgoing states are found to be eigenstates of the symmetry 
operator then they must be eigenstates with exactly the same eigenvalues, else this transition is not allowed. 

We begin the discussion with the parity operator V . We align the z axis along the direction of one of the incoming 
gluons. The outgoing or final state is parity-odd, as we know that our final state is the photon, or a state composed 
of a lepton and anti-lepton in the state. The gluons, on-shell in this calculation, are each parity-odd. We may still 
construct a parity-odd in state via the following method: we label the possible in states as 



|L+;L-), \L+-R-), \R+-L-). 

Where, the \R+-, R—) is the state where both gluons are right handed. The \L+;R—) state indicates that the 
gluon moving in the positive z direction is left-handed while that moving in the negative z direction is right handed 
(we have used the notation that the -I- sign indicates the gluon moving in the positive z direction). The parity 
operation interchanges the momenta of the two gluons but leaves the direction of their spins intact. Hence the state 
\R + R—) — \L + L—) is odd under parity operation. This implies that only this combination of incoming gluons is 
allowed by parity to fuse to form the virtual photon and hence the lepton pair. 

We next turn to the rotation operator, TZ. The in state is the state of two gluons; the out state may be considered 
to be either the temporary virtual photon, or the finally produced pair of lepton anti-lepton. One may chose either 
for this analysis; we decide on the photon as it is simpler. For the in state we use the only state that is allowed by 
parity i.e. \R-\- R—) — 1-^ + L—). This state may be re-expressed as the action of creation and annihilation operators 
on the vacuum state as, 



\R + R-) - \L + L-) 



4+4- 



|0). 



(38) 



Where |0) is the vacuum state. The creation operator a^_|_ creates a right handed gluon traveling in the positive z 
direction. The remaining creation operators have obvious meanings. The out state is the photon at rest and thus has 
the rotation properties of the spherical harmonics Yi.,„(0, 0). As the in state has both gluons either right handed or 
left handed, the z component of the net angular momentum is zero. Hence photon out state also has m = 0. 

We will rotate the in state and the out state by angle tt about the z-axis and then about the x-axis. The photon 
out state, mimicking the rotation properties of Yi,o(^, </>), is an eigenstate of either rotation with eigenvalues -1-1 and 
— 1 respectively. Focusing on the in state, we note that rotation by an angle (p about the axis n is achieved by the 
action of the appropriate operator U{R^) on the state in question. 



U{R%)\R+-R- 



UiR'l 
U{Rl. 



,t 



'R.+ \ 



ak+u-\R;)u{R;)a^^^u-\R;)u{R;m. 



(39) 



14 



Recalling the action of the rotation operators on creation operators (see Ref. [13), we obtain, 



(40) 



Where 'D{R^)fiii — {R\e^''''^\h) is the rotation matrix for the rotation of the state (in this case vector). The index 
h runs over all the possible z components of the spin of the particle. The vector p represents the new direction of 
motion of the particle after rotation. The action of any unitary operator, such as a rotation, on the vacuum will result 
in the vacuum again. Setting = tt we obtain the simple relation for the action of the rotation operator on the gluon 
creation operator. 



U{Rl)a\_U-\Rl) 



(41) 



Using the above it is not difficult to demonstrate that the in state of two gluons is an eigenstate of R^ with eigenvalue 
+1. Thus, both in state and out state are eigenstates of R^ with the same eigenvalue. As a result, there is no 
restriction to this transition, on the basis of this symmetry. 

We now concentrate on rotation by tt about the x axis. The outstate is an eigenstate of this operation with 
eigenvalue -1. Using Eq. H4U|I we note that, 

UiR:)a^j,.U-\R:) = 4_ 



UiR:)ai_U-\R:) = 
U{R^)al^U-\R^) = 
UiR:)a[^U'\R:) = 



''L+ 



(42) 



One may, thus, demonstrate that the two gluon in state is an eigenstate of the above rotation with eigenvalue +1, 

U{R:)i\R+;R~)~\L+;L-)) = [c/(i?:)al,+ C/-i(i?:)C/(i?:)at,„C/-i(i?^) 

- UiR:)a[^U-\R^^)U{R^M^U-\R;)] t/(i?:)|0) 

= \R-,R+) -\L-,L+). (43) 



This implies that this transition is not allowed by any interaction. Thus, we demonstrate Yang's theorem in the 
vacuum: this transition is not allowed 

B. Yang's theorem in media 

The above argument for no transition has been formulated for two massless vectors fusing to a spin one final state in 
the vacuum. We now intend to extend this to a transition in the medium. One may argue at this point that the correct 
method of analyzing this situation would be to start from a particular many body state; invoke the matrix element of 
the transition (this would give us the requisite creation and annihilation operators) and end up in a particular final 
many body state 



M 



f f f 



I / d'^xHiix) \n 



(44) 



This has to be followed by squaring the matrix element and weighting it by the Boltzmann factor e~^^', where Ei 
is the total energy of the in state, /3 is the inverse temperature. Then this quantity must be summed over all initial 
and final states to obtain the total transition probability per unit phase space for this process as 
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2 




(45) 



The above method though comprehensive, does not allow a simple amplitude analysis as the case for the vacuum. Such 
an analysis may be constructed by drawing on the spectator analysis of loop diagrams. This method as applicable to 
this process has been expounded in the previous section. Our results are essentially contained in Fig. 

The following analysis with spectators may appear to be rather heuristic at times. The reader not interested in 
such a discussion may consider the fact that the introduction of the medium formally involves the introduction of 
a new four-vector n into the problem. If we were to consider the case of dileptons produced back-to-back in the 
rest frame of the medium, the results from the vacuum should still hold as in this case the only new ingredient is a 
new four-vector of the bath (n = (1,0,0,0)). This four-vector is obviously rotationally invariant and cannot in any 
way introduce rotational non-invariance via dot or cross products with any three- vector in the problem. However, if 
the two gluons are not exactly back-to-back or equivalently the medium has a net three-momentum, then rotational 
invariance is explicitly broken. Even if we were to boost to the frame where the gluons are exactly back-to-back, we 
would find the medium streaming across the reaction. The above argument for the validity of the theorem for static 
dileptons will now be demonstrated via the spectator interpretation. 

We consider the Feynman diagrams of Figs. 01and|31 The effect of the medium, on the transition, is understood 
as a change in the in state to include an incoming quark from a particular quantum state a. Where, the index a will 
be used to indicate all the characteristics of the quark in question such as the momenta, spin or helicity, colour etc. 
The out state will also be modified as indicated to include a quark emanating from the transition and re-entering 
the medium in the same quantum state a vacated by the incoming extra quark. In the discussion that follows, we 
will keep referring to the original state containing the two incoming gluons as the in state, and the one with the 
outgoing dilepton as the out state. The extra particles that enter the reaction from the medium or exit the reaction 
and go back into the medium will be referred to as 'medium particles'. The full effect of the medium will only be 
incorporated on summation of the transition rates obtained by including all such states a weighting the entire process 
(incoming particles — > reaction — > outgoing particles) by appropriate thermal distribution factors for the incoming 
and outgoing medium particles. The thermal factors will essentially be those of Eq. H36|l . No doubt, there must also 
appear thermal factors for the incoming gluons. For the duration of the entire discussion, we will constrain the two 
gluons to have the same momenta; the distribution functions will thus play no role, and hence have been ignored. 

The new total in states and out states will now be given by state vectors that look like. 



In the above equation, we have taken the incoming and outgoing particle from the medium to be a fermion, as 
is appropriate in this case. The reader will recognize the thermal factors to be exactly those of Eq. (|36|) . Each 
state may once again be obtained by the action of the corresponding creation operators on the vacuum state. The 
new additional factors h{Ecr) are the appropriate distribution functions, used in the expressions to indicate particles 
leaving and entering the medium. Unlike the in and out states, the contributions from these medium states are added 
coherently, i.e. one does not square the amplitude and then sum over spins and momenta but rather the procedure is 
carried out in reverse, as indicated by Fig. |S1 The sum represents integration over all momenta, sum over spins 
and colours etc. 

Our method of extending Yang's symmetry will involve identifying certain subsets of the entire sum to be performed, 
which will turn out to be eigenstates of the rotation and parity operations to be carried out once more on these states. 
The argument will essentially be the following: if we can decompose the entire in and out state into certain subsets, 
with each subset being an eigcnstate of the symmetry operator with the same eigenvalue, then the entire in and 
out states will also be eigenstates with the same eigenvalues. Then, as for the vacuum process, we will compare the 
eigenvalues for the in state and outstate. 

To illustrate, we focus on a subset of four terms in the full sum in which one of the incoming mcdium-fcrmions has 
a three-momentum q. To keep the discussion simple we pick q to be in the yz plane (the discussion may be easily 
generalized to include g in an arbitrary direction). The four processes under consideration are: 



Y.[\R+; R-) \<J) -\L+;L-)\a)^ I7*) \<j) ) 1 2h{E,) . 



a 
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R+;R-)-\L+;L-))\q-A) 
(|i?+;i^-)-|L+;L-))|7^;rf;T) 
{\R+;R-)-\L+;L-))\n^q;i) 
R+,R-)-\L+;L-))\n%niq;i) 



|7*)k1T)J(l-2n(£;,-T)) 
|7*)|7^:q;T)](l-2n(£;7^.^,T)) 
\r)\niq-d)\{l-2fi{Enu,C)) 
\^*)\niKlq-i)\{l-2h{En^^nU-i))- 



(46) 



Where. TZ^q represents the three-momentum q rotated by an angle tt about the z axis, TZ^q represents q rotated about 
by an angle tt about the x axis. The arrows t, i represent the z component of the spin of the medium-fermion. As 
we are in the centre of mass of the thermal bath we have 



^qA — -^TC;«,T — ^n%q,i ~ ETiin%q,i- (47) 

Thus we may completely factor out the distribution functions. Without loss of generality we may combine all four in 
states and out states to give, 



(|i?+; R-) - \L+- L-)) \q- ]) + \nlq- ]) + \nlq- [) + IKTZlq: i) 
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|g; T) + mq; T> + \Kq; I) + IKUlq; I) 



(48) 



Now, it is simple to demonstrate using the methods of rotation of creation operators outlined in the vacuum case, 
that both the in and out states are eigenstates of 7?.J. Concentrating on the rotation of the in state we obtain 



[/(i?^) (|i?+; R-) - \L+; L-)) [ \q- ]) + \niq- ]) + {TZ^q; I) + miZ^q; I) 
U{Rl) (aJ,^+4^_ - 4;+4;-) U-\Rl) 

U{Rl)al,U-\Rl)U{Rl)al.^.^U-\R^^^^ \ |0) 

I (|i?+; R-) - |L+; L-)) \ |7^^<f; j) + |7^^7^^q; [) + |<f; T) + |7^;rf; T) 



(49) 



Note that the medium-fermions just mix into each other, but the over all state remains the same. Following the 
above method one can show that the outstate is also an eigenstate of but with an eigenvalue of i. Thus, we can 
decompose the entire sum over spins and integration over the three-momenta of the medium-fermions into sets of 
states as indicated, each will result in an in state and an out state between which no transition is allowed. For the 
rotation TZ^ we note that the eigenstates are in fact a subset of two states: in this case, the sum of the first two states 
of Eq. H46|) are eigenstates of 7?.^; as is the sum of the third and fourth state. 

This would imply that such a transition, as implied by the Feynman diagrams of Fig. 0] can not occur. There is 
however a caveat to the above discussion. Note that in the vacuum case we expressly boosted to the frame where the 
two gluons would be exactly back-to-back with their three-momenta equal and opposite. Then, rotational symmetry 
was invoked to demonstrate the impossibility of this transition. In the case of the processes occurring in medium, we 
tacitly began the analysis with the two gluons once again exactly back-to-back in the rest frame of the bath. However, 
if the two gluons are not exactly back-to-back or equivalently the medium has a net three-momentum, then rotational 
invariance is explicitly broken. Even if we were to boost to the frame where the gluons are exactly back-to-back, we 
would find the medium streaming across the reaction. This would make the distribution functions of the two gluons 
different (even though in this frame they have the same energy), Eq. H47|) would no longer hold. As a result it will 
not be possible to construct eigenstates of the rotation operators 7?.^ and TZ^^ as done previously. As the in and out 
states will no longer be eigenstates of TZ^ and TZ^ with different eigenvalues, transitions will, now, be allowed between 
them. 

In the above discussion, we have demonstrated how the medium may, once again, break another symmetry of the 
vacuum; in this case rotational symmetry. This allows the transition of Fig. |5] to take place in the medium. This 
process is strictly forbidden, in the vacuum, by two different symmetries (charge conjugation and rotation). It is 
forbidden in the exact back to back case by rotational symmetry in a C broken medium i.e. the effect is zero for 
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p = 0. To obtain a non-zero contribution, rotational symmetry has to be broken by a net p. The magnitude of the 
signal from such a symmetry breaking effect may only be deduced via detailed calculation. In the next section we 
shall outline just such a calculation. 

The derivation of Yang's theorem depended expressly on the the two incoming gluons being masslcss. This enforced 
their polarizations to be purely transverse. Another way of breaking rotational invariance is thus by giving masses to 
the gluons. This is not unjustified since in medium they acquire a thermal mass. If the gluons are considered as being 
massive, this implies that they now have three rather than two physical polarizations. The longitudinal polarization 
state is then physical and can be seen as being responsible for the breaking of Yang's theorem. Under parity, rotation 
around the x-axis of tt, and rotation around the z-axis of tt, the creation operator of the new 0-polarization state 
transforms, respectively, as 

Val^V-' = -al^ (50) 
UiR:)al^U-\R:) = -al^ (51) 
U{K)al±U~\K)^al^ (52) 

Then we can construct the two new in states with total angular momentum along the z-direction of +h: 

_L{|i?+;0-)-|0+;L-)} (53) 

and —h: 

-^{|0+;i?-)-|L+;0-)} (54) 

We must now inquire as to the possibility of a transition with a virtual final photon in the m = ±1 states. First note 
that now the out states are not eigenstates of U{R%) operator, but still are of U{R^) with eigenvalue —1. Applying 
this operator on the in states, one finds that they are eigenvectors with eigenvalue —1. Therefore, the in states and 
out states share the same eigenvalues. Thus the transition is not prohibited by parity and rotational symmetries. 

In a real medium both effects discussed {i.e. finite momentum and massive gluons) would be present and simulta- 
neously lead to the breaking of this symmetry. In this article we separate these two effects and study each in turn. 
The complete calculation incorporating the two simultaneously will be left for a future effort . 

There remains yet another means by which the symmetry of Yang's theorem may be broken: that of a rotationally 
non-invariant regulator. The results from this scenario have already been presented in Rcfs. [^Q- In these calculations 
the dileptons were produced back-to-back, i.e. from a virtual photon which is static in the rest frame of the plasma 
from the fusion of massless gluons. The reader will note that Yang's theorem holds for such a process and hence should 
result in a vanishing rate. Yet this rate was found to be non-zero. The reason behind this result is the choice of the 
regulator used in those calculations. From Eqs. H48() and H49() we note that we required the coherent sum of at least 
four quark states, with the incoming quark occupying symmetric angles, for Yang's theorem to hold. If we designate 
one of the gluons to be along the z axis, and one of the incoming quarks is assigned the momenta {q,6,(l)), then the 
configuration that obeys the rotational symmetry of Yang's theorem will include incoming quarks at {q,9,(j) + tt), 
{q, TT ~ 9,(p) and (g, tt — 0, + tt). Thus in the 9 integration, one must include balancing contributions 9 and tt — 9. 

As is the case in this article, the results of Refs. 0:13 consisted of the sum of contributions from multiple residues, 
some of which displayed singularities as ^ or ^ tt. This corresponds to the one of the internal lines in Fig. Agoing 
on shell. This divergence is canceled when all the different residues are combined, as will be shown in the next section. 
Each residue is then evaluated using a regulator. Two obvious choices are the angle S < 9 < t: ^ 6, and the magnitude 
of three-momentum of the intermediate state x = \q + k\ ^ q'^ + k'^ + 2kq cos 9 where \q — k\ + e < x < q + k — e. 
All residues are then evaluated in the limit 5 ^ or e ^ where the same regulator is used throughout. The 
divergence will be canceled in either case when all the residues are summed. In the calculations of Refs. 0, Q the 
three-momentum regulator was chosen. 

From the preceding discussion, it is obvious that integration using the angular regulator S obeys Yang's theorem as 
symmetric contributions from 9^6 and 9 = n—d are included. However this is not the case with the three-momentum 
regulator. Though e = corresponds to 5 = 0, at e ^ we find that x ~ \q + k\ ~ e corresponds to a, 9 — 5i ^ 0, 
while X = \q — k\ + e corresponds to a = tt — ^ 0. After some calculation, it may be demonstrated that 
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while 



Thus for a given k,q the same e corresponds to different hmits of the angular integration. The rotational invariance 
required for Yang's theorem is broken and a non-vanishing contribution results. The physical interpretation of the 
non- vanishing results of Refs. Q thus become unclear. The results obtained here do not depend on those prior 
findings. In this sense, the current article constitutes an update and a correction. In what follows, the symmetry in 
Yang's theorem will be broken only by real physical effects prevalent in hot media. 



V. RESULTS FOR p / 

We concentrate on the breaking of the symmetry underwriting Yang's theorem by imposing a non-zero three- 
momentum to the process. In other words, the virtual photon has a net three-momentum in the rest frame of the 
bath. The gluons are considered to be massless. As may be easily understood from the preceding section, the 
magnitude of symmetry breaking rises with the magnitude of the three-momentum. Thus the largest possible values 
of p will lead to the largest signals. 

We want in the end to calculate the number of dilcptons per unit spacetime per unit energy per unit \p\ ~ p, i.e., 

dxHEdp ~ dEdp ~ j It^ 

Where we have integrated the differential rate over all solid angles fij,. The angle of p is always measured from the 
direction of the more energetic incoming parton: in this case the gluon with energy > iJ/2. This procedure will also 
be followed for the Born term and will be explained in greater detail in the last subsection. 

We begin by presenting results for a simpler case. We look at the differential rate when the two incoming quarks 
or gluons are forced to be back-to-back but have different energies. This may be obtained by setting the angle 5 
between the photon and the in coming gluon to zero. Alternatively one may obtain this by expanding the rate in a 
Taylor expansion in angle and keeping only the first term. One reason for considering this special case is that this is 
the simplest generalization from the p — Q case. Yet another reason for considering this case is the possibility of an 
analytical solution. We provide complete analytic results in this case, as opposed to the general case where the final 
integrations can only be performed numerically. 



A. Differential rate for back-to-back gluons 

We begin by evaluating the sum of the six matrix elements on the right hand side of Eq. H37|) . Recall that i,j 
indicate the polarizations of the incoming gluons where as I is the polarization of the outgoing photon. In this 
configuration, a variety of simplifying relations result: 

M+,-,i=M-^+,i = G (57) 

as expected (and pointed out before), in a back-to-back situation, both gluons have to arrive with the same polarization, 
i.e. either both must be right handed or both left handed. This will result in a net z component of angular momentum 
Lz ~ 0. The configuration with one left-handed and one right-handed in a back-to-back scenario will have an = 2 
and thus will not couple to a spin one object. 

= 7V/(,j.- = (58) 

There is no contribution to the transverse modes of the photon. This is obvious as the two gluons form a system with 
a net = 0. 

We are now in a position to write down explicit expressions for the various matrix elements of Eq. (|37|) . i.e. Mi.j^i ~ 

+ + + t^i,3,i + f^i,j,i t^i,3,i) ^'^(-^9' A*)- performing aU the angular integrations, frequency sums 
and contractions with polarization vectors, the results are: 
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(59) 

We are thus concentrating on the virtual photon produced only by back-to-back gluons of unequal momenta, and 
will compare with the rate of production from only back-to-back quarks of unequal momenta. We are thus breaking 
Yang's symmetry by the introduction of a net three-momentum p. It is a simple exercise to check that the above 
matrix clement vanishes linearly with p as p ^ 0. The apparent pole in g, is canceled between the six terms. There 
is still the dq integration to be performed, this is done numerically. 

The differential production rate for pairs of massless leptons with total energy E and and total momentum p is 
given in terms of the discontinuity in the photon self-energy as (see Eq. 



dR _ r 1 

~ 3(27r)5 (p2) _ l 



(60) 



Where r 



Disc [— in^'']. In general, the matrix clement depends on the angle 5 between p and k. As a 



result, two equivalent means of angular integration present themselves: we may set p along the z direction; in which 
case the matrix element will depend on the polar angle 5 = Ok oi k. As a result, the integral over Op yields an overall 
factor of 27r. Alternatively we may set k along the z direction, in which case the rate depends on the angle 5 = Op, 
and the integral over 51^ yields an overall factor of 2-k. Both methods are equivalent. We chose the latter prescription. 
Thus we calculate the derivative of the differential rate with respect to the incoming gluon angle ^Ik 



dR , . 

d'^pdnk 3(27r)5(p2)[e'9^^-l]df7fe' ^ ' 

As mentioned beforCjtemperatures in the plasma formed at RHIC and LHC have been predicted to lie in the range 
from 300-800 McV [HEl- For this calculation, we use T ^ 400 McV and 800 MeV. To evaluate the effect of a finite 
chemical potential we perform the calculation with two values of chemical potential fi — O.IT (left plot in Fig. ISJ and 
/i = 0.5T (Right plot in Fig. ^ |24l |. This calculation is performed for two flavours of quarks with current masses. 

In Fig. El the differential rate (Eq. for the production of dilcptons with an invariant mass from to 156 McV 

is presented. The energy is held fixed at 500 MeV and the three-momentum p of the dilepton is varied (a dilepton 
invariant mass M = 156 MeV for an energy E ~ 500 MeV corresponds to the three-momentum of the dilepton 
p = 475 MeV). In the figures, the dashed line is the rate from tree level qq ; the solid line is that from the process 
gg e+e~. We note that in both cases the gluon-gluon process dominates at very low mass and dies out at higher 
mass leaving the qq process dominant at higher mass or lower momentum. The Born term displays a sharp cutoff 
at photon invariant mass M = \j2mE. The back-to-back annihilation of two massive quarks (of mass m) to form 
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FIG. 6: The differential production rate of low mass dileptons from two back-to-back processes. Invariant mass runs from 
156 MeV to MeV. The energy of the dilepton is E — 500MeV, and the abscissa is the three-momentum p. The dashed line 
represents the contribution from the process qq e^e~ . The solid line corresponds to the process gg e^e~ . Temperature 
is 400 MeV. Quark chemical potential is O.IT. The second figure is the same as the first but with fj, = 0.5T 



the virtual photon of energy E and invariant mass M is no longer kincmatically allowed. Also, the annihilation of a 
quark antiquark pair to form a dilepton is not allowed for any incoming angle for dileptons with an invariant mass 
M < 2m. The gluons being massless, continue to contribute in this region: this contribution is shown in the right 
panel of Fig. [7| Thus the signal from gg fusion (for partons with current masses) is dominant at low invariant masses 
for intermediate dilepton energies. In Fig. we indicate the influence of a higher plasma temperature on the rates. 
Here, a plasma temperature of 800 MeV and /i = Q.hT is used; the left panel displays the rates below the Born term 
threshold and the right panel displays the rates above threshold. We note in the left panel of Fig. [7| as expected, 
that the giuon fusion term rises further due to thermal loop enhancement, however this rise is rather minimal. In the 
right panel we note that the rates for gg e+e^ continue to rise due to the growing distribution functions for soft 
gluons. In a realistic calculation the gluons would be endowed with a thermal mass that would cut off the steep rise 
in the gg rate. 
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FIG. 7: Left panel is same as Fig. |H|but with a temperature of 800 MeV. Quark chemical potential is 0.5T. Right panel is the 
rate of gg — > e'^e~ just below and beyond the Born term threshold. Born term threshold is at p = 499.6 MeV for a dilepton 
energy of _E = 500 MeV and quark current masses of 10 MeV. 
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B. Full differential rates 



In the previous subsection, we calculated the six matrix elements of Eq. (|37|l (or Fig. 0J, for the special case of the 
angle between k and p i.e. ^ = 0. A close inspection of the diagrams of Fig. 0]will convince the reader that, in the 
event that the in coming gluons and the outgoing photon are in the same line, the angular integration over all quark 
directions displays an azimuthal symmetry. This implies that the integration d(j)g results in a mere overall factor of 
2tt. The d6 integration, though non-trivial, can be performed analytically and results in Eq. H59|l . The integration 
over the quark magnitude cannot be performed analytically and we resorted to numerical means. 

In the case of a (5 ^ 0, the azimuthal symmetry in the cPq integration is absent and both d(j)q and d0q integrations 
are non trivial. It is no longer possible to perform both analytically. Following the dq integration, we also have 
to perform the integration over the angle S. This will give us the differential rate -jgj^. This turns out to be a 
complicated problem to solve in general. However, from the previous subsection we have learnt that the rate form 
this process is comparable to the Born term only at very low mass M = ^ — j? — > 0, or rather p ^ E (see Figs. 
1617(1 . If we insist on calculating solely in this limit an approximation scheme may be constructed. 

There are two basic scales that we input into this problem: the mass of the quarks m and the temperature of the 
plasma T (the chemical potential is always estimated as a fraction of the temperature: hence it does not constitute a 
separate scale). At the low invariant masses (of observed dileptons) in question the strange quark does not contribute. 
For the up and the down quark, we are considering plasmas where m << T. We now insist on observing dileptons 
with large four- momenta E,p ^ T, yet very small invariant mass AI = \/ E'^ — p^ ^ m << T. Yet another smaller 
scale is that oi x ^ E — p where a; = ^ ^ « M. One may construct three dimcnsionless scales from these 
quantities: 1 >> ^ >> ■§. 

We denote the incoming gluons by their polarizations Now, say i is more energetic and is ascribed the momenta 
k, the other (j) has momentum \p — k\ by conservation. As outlined in the previous section, we intend to integrate 
5 from (5 = when the two gluons are back-to-back and k » \p—k\ = E — k (lisa very hard gluon and j is very 
soft) up to (5 = Smax whcrc k = |p — k\ ~ E/2 (where, throughout the gluon i is more energetic than the gluon j). 
The remainder of the S integration may be obtained by simply replacing i with j and noting that the remainder is 
nothing but the same integration with the gluon j now ascribed the larger energy k and 5 defined such that the j 
points in the positive z direction. The magnitude of 5max may be estimated simply from the preceding discussion. 
As the gluons are massless 



E-k^ E^P_j:^ = \p- k\ = + fc2 _ 2pk cos{5) 



(62) 



This implies that 



^Ek - 

cos{S) = — (63) 



The value of (5 = Smax occurs a.t k — E/2, hence 



C0s{6,nax) = 1 

= 1 

Thus S^ax = - §■ 

Yet another inference about the behaviour of the rate with 5 may be drawn from Eq. H32|l . Here we note the 
presence of two Bose-Einstein distribution functions: n{k),n{E — k) and the factor in the measure. Note that as 
(5 — > 0, fc tends to its maximum value, and tu ^ E — k tends to its minimum value. This greatly enhances the factor 
J'_y, n{k)n{E — k), as compared to its value at S = Smax, where k and E — k are of the same magnitude. This implies 
that if the the matrix element does not rise sharply with 5 then the differential rate falls off as 5 is raised. 

The above mentioned observations allow us to expand the matrix element in a series in 5. On expansion we note 
the following behaviour for the longitudinal photon : 



6 p 
Y ^ E 

X 

E 



(64) 



Mi^j^siS) = mo - 7712(5^ + 
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Where the m^'s are all positive contributions, and depend on E,p,T,ii. Thus for small 6 the matrix element (or the 
square of the matrix element) drops as S rises from zero. A plot of this behaviour (in arbitrary units) for a typical 
case is shown in Fig. |H1 It should be pointed out that the angular integrations of the quark momenta {i.e. 6,(j)) 
may be analytically performed only after the expansion in S. The remaining integration over the quark momentum q 
is performed numerically. If the matrix element were not expanded in a series in S one would have to perform four 
sets of integrations numerically. The presence of poles in the matrix elements of the diagrams in Fig. 01 makes this a 
prohibitively difficult procedure. 




0.003 




0.003 



FIG. 8: Left panel shows the behaviour of (M = mo — m2S^)/M{S = 0) as a function of S for a typical case of T = O.SGeV, 
fi = 0.5T, E = IGeV and p = 0.9999GeV. Right panel shows the behaviour of Im(n'^)/Im(n'^){S = 0) as a function of S. 



Now that we have moved away form the back-to-back scenario, we will also witness the production of transverse 
photons. On expansion in S we note the following behaviour: 



Mij^±{S) = mi5 - 77138^ + ... 

As expected this contribution goes to zero as J ^ 0. It also turns out to very much smaller than the longitudinal 
contribution in the limit of small invariant mass dileptons. We thus obtain that the dominant contribution to the 
rate emanates from the longitudinal photons. 

It should be pointed out that such an expansion is only strictly valid as long as S is of the order of the smallest 
scale in the problem i.e. ^. However, as demonstrated by the plot of the imaginary part of the self energy in Fig. 
|51 it remains valid much beyond this point. This is due to the influence of the measure and the gluon distribution 
functions which drop rapidly as one moves away from S = Q. It should be pointed out that only the matrix elements 
have been expanded in a series, all other factors (e.g.ghion distribution functions, measures etc. ) retain their closed 
expressions in S. 

We are now in a position to integrate over 6. This is also performed numerically. As we have expanded the matrix 
elements in S and retained only a finite number of terms, the square of the matrix elements \A4\'^ grow beyond a 
certain S = 61. This growth is not real and is merely a facet of our finite expansion in i5. Two possible means of 
carrying out the S integration present themselves. We may terminate the integral at i5i. Ostensibly, this represents 
the lower limit of the rate. These are represented in Fig. |51as the solid lines (both with and without symbols). We 
may continue to integrate up to S = Smax] this will include the integration of a growing rate convoluted with an 
increasing angular measure (sin (5). This represents the upper limit of the rate. These are represented in Fig. |51as the 
dotted lines. As the invariant mass is lowered or the energy of the pair is raised the differential rate drops sharply 
from its value at 6 = 0; this invariably results in 7712 >> mi. As a result the integration beyond S ~ Si produces a 
large contribution; this in turn leads to the excess growth displayed by the upper limit at small invariant mass and 
large energies. 

The rates after integration over 6 are presented in Fig. |5| The dot-dashed lines arc the rates from the Born term 
at various energies of the virtual photon. The solid lines are the lower limits of the rates from gluon gluon fusion for 
the respective energies of the virtual photon. The dotted lines are the upper limits of the respective rates. In general 
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FIG. 9: (Color online) Differential rate for the production of dileptons at high momenta and small invariant mass. See text for 
details. 



the rates from gluon gluon fusion are suppressed as compared to the Born term except at very high momenta or very 
low invariant mass. Three cases have been presented where the energy of the dilepton is set at 0.25, 0.5 and IGeV. 
As may be noted, the results are quite similar to the back-to-back gluon fusion rates, i.e. the case at (5 = 0. The 
rate from gluon fusion rises beyond the Born threshold due to the rising Bose-Einstein distributions of soft massless 
gluons. 



VI. RESULTS FOR mg > 

We now turn to the case where gluons acquire a medium-induced mass and where the virtual photon has no net 
three-momentum {p — 0). As we have seen, the gluon's longitudinal degree of polarization allows to circumvent Yang's 
theorem. More specifically, wc expect transitions to occur when the net component of angular momentum is ±h: 



M,j,3 = 
M±,3,± = A^3,±,± = 



Trig eg 5 



2sbc 



dq 



■An{Eg,^i)J{q,k) 



(65) 
(66) 
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where 



where = \J k'^ + t^^^g- Therefore, in the back-to-back configuration with massive and equaUy energetic gluons we 

find that the non-zero matrix elements are proportional to the gluon mass and scale like ~ fi confirming the finite 
density nature of this process. We point out that the integrand diverges unless < 2mq. Beyond this limit the 
self-energy analytical structure becomes intricate. In this section, we present results only where the above relation 
holds. Dileptons production rates from regions beyond this threshold, as well as rates emanating from the fusion of 
massive gluons with p 7^ will be addressed in a future effort . 

To explore dilepton production in the range where the mass inequality is respected, we begin by ascribing current 
masses of 10 MeV to the quarks, setting the quark chemical potential to /i = 40 MeV, and the gluon mass to almost 
twice that of the quark (i.e. nig — 19.99 MeV) as a reference point. With these, we see that the differential production 
rate due to the gluon fusion is lower than the contribution from the Born term across the range of invariant mass 
(see Fig. llO|l . However, if we increase the quark chemical potential to /i = 200 MeV or reduce the quark current mass 
to niq = 1 MeV while maintaining the gluon to quark mass ratio {i.e. 1.999), we see that the gluon fusion rate may 
dominate over the Born term up to an invariant mass of 125 MeV in the case where the quark chemical potential 
reaches 200 MeV. We also present results where the gluon mass is set equal to that of the quark (i.e. 10 MeV), in this 
case we see that the rate from the gluon fusion dominates at low invariant mass. If we lower the gluon mass beyond 
that of the quark then the process will have a non- vanishing contribution in a region forbidden for the qq process due 
to its threshold (see Fig. [TT|l . 

If instead of current masses, the quark masses are set to values of the order of gT then we find the rate from gluon 
fusion to be suppressed as compared to the Born term (see Fig. lllH . This is not unexpected as the gluons fuse through 
a quark triangle; the presence of large masses in the quark propagators leads to the rates being suppressed in this 
region of parameter space. 

As in the case of gluon fusion with p 7^ 0, an interesting feature of the differential rate is its strong dependence on 
jj? (Fig. I12|) and its weak dependence on temperature and other energy scales: as the chemical potential increases, 
the rate rises. For the Born term the opposite behaviour is true. As the chemical potential increases, the antiquark 
population is depleted, inhibiting the production of dileptons through this channel. Thus, an accurate estimate of the 
differential rate will require a good knowledge of the baryon chemical potential as well as its variation with time in a 
QGP. 

VII. DISCUSSIONS AND CONCLUSIONS 

In this article we have presented a detailed study of the observational effects of broken charge conjugation, and 
broken rotational invariance in a QGP formed in a heavy-ion collision. The signal under consideration was the 
spectrum of dileptons emanating from such a medium. The reason behind this choice is evident: electromagnetic 
signatures provide a direct probe of all time sectors of a heavy-ion collision. The breaking of these symmetries is 
manifested at lowest order in the spectrum of dileptons produced by two gluon fusion into a virtual photon through 
a quark loop. 

Such a process is forbidden in the vacuum by both Furry's theorem and Yang's theorem. Charge conjugation was 
broken explicitly by the introduction of a non-zero population of u and d valence quarks in the plasma. A non-zero 
baryon density, present solely in these flavours causes a net electric charge density in the medium. This leads to the 
breaking of Furry's theorem which holds purely in neutral media. 

The presence of a preferred rest frame of the medium leads to the introduction of a bath four-vector n into the 
problem. If calculations are performed in this frame then n= (1, 0, 0, 0). No such vector exists in the vacuum. If two 
massless back-to-back gluons with equal energy fuse in vacuum along the z-axis through the quark triangle, then the 
out state consisting of a static virtual photon will have the z component of its spin = 0. Both in state and out 
state are eigenstates of a rotation by tt about the a;- axis but have different eigenvalues. Hence such a transition is 
forbidden, this is the statement of Yang's theorem and is based on the invariance of both in state and out state under 
a rotation by tt about the x-axis. 

The above argument of no transition (due to the in state and out state possessing different eigenvalues with respect 
to the rotation by tt about the x-axis) continues to hold even for the production of a static photon in the rest frame 
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FIG. 10: (Color online) DiflFerential production rates for T = 400 MeV. See text for details 



of the bath where the in state and out state are modified to include spin-half fermions. It should be pointed out 
that we have tacitly assumed the plasma to be infinite in extent and isotropic; a realistic plasma of finite extent 
which is not spherically symmetric will explicitly break rotational invariance. To our knowledge, this fact may be 
understood solely in the spectator interpretation of loop diagrams. The spectator interpretation represents a formal 
procedure by which the imaginary part of a diagram containing loops may be re-expressed in terms of the product 
of matrix elements consisting solely of tree diagrams and particles from the bath that do not partake in the reaction 
process. The spectator interpretation for the imaginary part of a three-loop diagram was derived for this process and 
is essentially contained in Fig. |S1 In the spectator interpretation different states containing fermions with different 
spins are added coherently. This allowed us to construct a subset of the entire sum of in states that respected the 
rotation symmetry of the vacuum state (see Eq. H48|l '). For each such subset no transition was allowed by arguments 
similar to those used in the construction of Yang's theorem (see Sec. IV B). 

This invariance will be broken by the presence of any three vector in the problem. The two possible choices for such 
a three- vector are a net three- momentum of the virtual photon {p ^ in the bath frame), or a non-zero z-component 
of its spin {Jz ^0). In the first case one may boost to the rest frame of the static photon. However, this will no longer 
be the rest frame of the plasma and rotational invariance will be explicitly broken. In the second case the production 
of a virtual photon with a ^ will require an incoming massive gluon which breaks one of the principal conditions 
required for Yang's theorem to hold. In a real QGP we would expect both effects to be present simultaneously. In the 
interest of a clearer understanding of the mechanism of symmetry breaking we had chosen to explore both possibilities 
in isolation. 

In the case of a virtual photon with a net three-momentum p, we began with the case of two massless gluons in a 
back-to-back configuration along the z-axis but with unequal energies. This resulted in a virtual photon with a net 
three-momentum along the z-axis. Under this kinematic restriction we present a closed analytical expression for the 
matrix element (see Eq. (jSHJ)- The rates from this matrix element are plotted in Figs. El and din comparison with 
the Born term. These figures represent the cases with plasma temperatures at 400 MeV and 800 MeV, corresponding 
to the cases of QGP formed at RHIC and LHC energies. We find as expected that the differential rate rises with 



26 



10 



10 



10 

10-' 

^ , 

10-' 



10 



10 



500 



Bom Term 




gg->e+e-: m^= 


1.5m 




q 


.-- gg->e+e-:m^ = 


0.5m 







1000 



1500 

E (MeV) 



2000 



2500 



FIG. 11: Differential production rates for 2g ^ II with T = 400 MeV, fiq = 40 MeV, m, = 250 MeV, and different gluon 
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increasing chemical potential /i; two extreme cases with /i = O.IT and fi = 0.5T have been explored. The rates however 
show a rather modest rise with increasing temperature. This has been pointed out earlier due to the vanishing of the 
g^T^ component in the HTL calculation of this loop diagram 0. The rates are also observed to rise with increasing 
three-momentum \p\ (for a given fixed energy) as expected due to the breaking of Yang's theorem. As the gluons 
are massless they continue to contribute into the region beyond the Born threshold (see Fig. |7|). In this region, 
contributions originate from the fusion of a gluon carrying a large majority of the photons energy with an ultra-soft 
gluon carrying a tiny fraction of the photon energy. In fact the rising rate in the right panel of Fig. is due to 
the Bose enhancement obtained from the distribution function of the soft gluon. The presence of a gluon dispersion 
relation or a gluon mass will lead to this rate reaching a maximum at a threshold set by twice the gluon mass. 

Integrating over all incoming gluon angles turned out to be an involved procedure and led us to invoke an approx- 
imation scheme. At very small invariant mass, we noted that the gluon fusion rate is dominated by back-to-back 
gluon fusion. We thus expanded in the angle between the photon and the more energetic gluon S. Results for the 
differential rate per unit energy and per unit momentum d?R/dEdp have been presented in Fig. |51 in comparison 
with the Born term. As expected the rate from gluon-gluon fusion becomes comparable to the Born term only at 
very low invariant mass, for photon energies of 0.25, 0.5 and IGeV. As noted in the previous section a large portion 
of the enhancement may be attributed solely to the lack of a mass for the gluons and Bosc-Einstein distributions as 
opposed to a Fermi-Dirac distributions for the quarks. 

In the case of a virtual photon with a net J^, the possible choices are Jz = ±1. This requires one of the incoming 
gluons to be in a longitudinally polarized state. Hence, the gluons were endowed with a mass. One may ascribe the 
origin of such a mass to dispersion in the medium. As in the previous case, the rate is seen to rise sharply with 
increasing chemical potential. Due to analytic considerations, the quark mass (rUg) was always set to be larger than 
half the gluon mass (nig < 2mg). In this kinematic region, the rates from gluon gluon fusion turned out to dominate 
over the Born term for low invariant masses of dilcptons, if the quarks were chosen to be light nig << T. However, 
the rates were sub-dominant to the Born term for the production of dileptons with large invariant mass, or for quark 
masses ^ gT. We point out, that, in this calculation T = AQQMeV , hence, g ^ 2. Thus unlike the case for plasmas 
at very high temperature gT ~ T. 

Gluons with masses at and above this threshold (2mg) may decay into two quarks which arc both simultaneously 
on shell. In the language of spectators, this corresponds to one of the propagators in the diagrams of Fig. 01 going on 
shell. The calculational and interpretive complications that arise from this situation are rather involved and represent 
a problem for the spectator interpretation. A preliminary calculation without the use of the spectator interpretation 
at this threshold in the limit of massless quarks and gluons found the rate to be large ; however these required the 
use of momentum cutoffs which did not respect the symmetry required by Yang's theorem. As pointed out earlier 
this made the physical interpretation of these results unclear. The computation of the rates at and beyond this 
threshold will full quark and gluon dispersion relations at finite three-momentum and their interpretation in terms 
of the spectator picture will be dealt with in a subsequent calculation |2^. Our goals in the present article have 
been to separately elucidate certain symmetries of the vacuum which are broken by a particular channel of dilepton 
production in a quark gluon plasma. For the purposes of simplicity we computed the rates from this channel for a 
plasma in complete thermal and chemical equilibrium. As the dileptons in this channel are produced essentially form 
the fusion of gluons, this process may display far greater significance in early plasmas which are estimated to be out of 
chemical equilibrium with large gluon populations. An accurate extimation of the full dilepton rate from this channel 
will require a two step process. One needs to combine the two effects of symmetry breaking discussed in this article. 
This rate will then have to be folded in with a realistic space-time model of the evolution of the plasma. Such a model 
will have to include an estimation of the early gluon population with estimates for the effective inmedium masses of 
the gluons and the evolution of these quantities with time. Work in this direction is currently in progress [26|. 
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APPENDIX A: CONTOUR INTEGRATION OF T""^" 

111 this appendix, wc outline the formal calculation of the two gluon photon vertex in the imaginary time formalism, 
using the method of contour integration. In this case, the standard method of contour integration will be modified to 
allow for the appearance of expressions which may be easily generalized from the case at zero density. This procedure 
allows for the construction of the spectator interpretation. As mentioned before in Eqs. Q and (O the Feynman 
rules for the two-gluon-photon vertex in the imaginary time formalism are. 
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Where the trace is implied over both colour and spin indices. The zcroth components of each four- momentum are 



(A2) 



q° = i(2n + l)7rT + /i. 



i2rm:T, 



fc" = i2jTTT. 



where n, m, j are integers, fi is the quark chemical potential. The overall minus sign is due to the fermion loop. The 
sum over n runs over all integers from — cxj to +00. As mentioned previously, the momentum dependent and mass 
dependent parts of the numerators of Eqs. HAlj) and HA2|) may be separated as 
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Where y^^*''''''' represents the trace of four 7 matrices and S^'"'''3p7 represents the trace of six 7 matrices. Employing 
the methods of residue calculus, the sum over n may be formally rewritten as a contour integration over the infinite 
set of contours Ci (See Fig. [T^ 



TE/(«" = *(2n+l)7rr + M) = ^ / dg"/(9°)^/3tanh Q/?((zO - a^)) (A5) 

The contours Ci may be deformed to those of C2 (see Fig. These are a set of two linear contours meeting at 
±100, one from q^ = —ioo+^+e — > g*^ = ioo+^+e, and another from ioo+^—e — s- —ioo+^—e. Here, and henceforth in 



29 







i ♦ 




QJ Tt i(3)T+ 11 


• Tt i(3)T+ [i 




C^) }■ TT in \T + u. 


■ , JT i(l)T+ ^ 




Qj ni(-l)T+ ^ 


' . ni(-l)T+ ^ 




QJ Ki(-3)T+ H 


. ni(-3)T+ (i 



Contour C i 



Contour C 2 



tKi(l)T+ ii 



ni(-l)T+ P-\ 
ni(-3)T+ H \ 



Contour C 



Contour C 



3 ^^^^^^uK ^ 4 

FIG. 13: The contours used to evaluate the Matsubara sum with a finite chemical potential. 



See text for details. 



all discussions of contours, residues and analytic continuations, e will represent a vanishingly small quantity. One may 
now proceed by the standard method of and separate a vacuum part, thermal part and a pure density contribution. 
Instead, another set of contours is introduced: these run along the y axis from = —ioo + e — > g*^ = ioo + e and from 
ioo — e — > —ioo — e. Admittedly, as e ^ this contour will produce a vanishing contribution. The integrand in Eq. 
(|A1|) has six powers of in the denominator and only three in the numerator. Hence it vanishes faster than a linear 
term as 00. As a result, this quantity obeys Jordan's Lemma and the two integration contours around and /i 

may be connected by line segments at ±ioo. These line segments shown as curved lines in the third contour of Fig. 
^|will have zero contribution to the entire integral. The total contour thus obtained is referred to as C3. We now 
split the integrand into two, one piece for all the contours on the positive of the x axis denoted as C3 , one piece for 
the sole contour on the negative side of the x axis denoted as C|, i.e., 



I- £^ .."/(.")i/3tanh - M)) = ^ J~^~[ ^."/(.^ ^ ' 

Co 



- ^r.„o_,, , . ■ (A6) 



27ri Jioo+^i-e J -icxi+fj,+e/ Qa \2 e^*'? -|- 1 

The terms may now be separated into a vacuum piece and a matter piece, note the similarity between this and the 
zero density separation. In this procedure, we differ from the standard method Q in not extracting an explicit finite 
density piece. The main reason for the extra contour deformation is to obtain the final answer in a form where the 
zero density contribution is obvious. In this spirit, we now reverse the direction of integration in C| and note that 
the vacuum piece has no poles at i{2n + l)7rT + fi. Thus the contours in the vacuum term may be allowed to overlap 
by setting e = 0. We obtain 
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We now let e — > on the contours on the positive side of the x axis. This procedure will deform the two linear contours 
at /X ± e back to the small circles around the points i{2n + l)i:T , this part will become similar to the initial contour 
Ci. The rest of the contour can be closed by including the infinite arc in the = +oo direction in the clockwise 
sense. This multiply connected contour is indicated as Cf — Ci and displayed on the right of the fourth plot in Fig. 
1131 The linear contour on the negative side may be closed off as always by the infinite arc extending to q^ = — oo. 
This is indicated as C\ and shown as the left contour in the fourth plot of the figure. The contour integration over 
either contour may be replaced by the sum over all the residues at all the poles enclosed by the contour. Note that 
the poles at i{2n + Vji^T + /i, excluded by the multiply connected contour, arc not to be included in the sum over 
residues. Thus our final, formal result is. 
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We may substitute the full integrand in Eq. (|A1(I to obtain the result of contour integration as 
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A similar contour analysis as above may be performed for Eq. (|A2|I . with the added extra step of setting q*^ — > 
— g*',(f —q. This procedure will produce a final contour of integration which is a mirror image of C4. There will, 
once again, be an infinite semi-circle extending to +00 connected with the line running from —ioo + e ioo + e. 
There will also be an infinite semi-circle extending to —00 connected to the vertical line running on the negative side 
of the X axis. This contour will however be multiply connected with the poles at —i{2n + l)7rT — /x excluded from 
the region bounded by the infinite semi-circle. As before these poles shall be excluded from the sum over residues. 
Following this procedure, one obtains the result of the contour integration for Eq. (|A2p as. 
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Note, that the vacuum term is at least naively linearly divergent and thus the shift in momentum integrations may 
not be performed as above. However, ft'om Furry's theorem, one obtains that the sum of the vacuum terms from 
Eqs. HA9|I and (| A 10(1 must be identically zero. Also note that the presence of the thermal distribution functions over 
quark momenta renders these integrals ultra-violet finite. Quark momentum shifts are thus definitely allowed for the 
thermal parts of Eqs. (|A9|) and IjAlOp . Hence, We ignore the vacuum pieces and combine the matter pieces of both 
terms to obtain T^"'P = T^"'P + Tf^P" as. 
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